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O ABSTRACT 

We report the detection of Lyman-Werner absorption from molecular hydrogen (H2) 
I ^ I at z — 0.56 in a sub-damped Lya system with neutral hydrogen column density A'hi ~ 

2^q19.5±o.2 cjjj-z rpjjjg j.jjg ^j-gj. system analysed at a redshift < 1.5 beyond the MiUcy 
I Way halo. It has a surprisingly high molecular fraction: logio /h2 > —1-93 ± 0.36 based on 

J> modelling the line profiles, with a robust model-independent lower limit of > 10~^. This 

•/^ is higher than values seen along sightlines with similar iVni through the Milky Way disk, 

the Magellanic clouds, or towards most higher redshift QSOs. The metallicity of the absorber 
is 0.14_Q jg solar, with a dust-to-gas ratio 0.041q q3 of the value in the solar neighbourhood. 
Absorption from associated low-ionisation metal transitions such as O I and Fe 11 is observed 
in addition to O VI. Using CLOUDY models we show that there are three phases present; a 
^-H ^ 100 K phase giving rise to H2, a ^ 10* K phase where most of the low-ionisation metal 

absorption is produced; and a hotter phase associated with O VI. Based on similarities to high 
velocity clouds in the Milky Way halo showing H2 and the presence of two nearby galaxy can- 
J> didates with impact parameters of ~ lOkpc, we suggest that the absorber may be produced 

" '"j by a tidally-stripped structure similar to the Magellanic Stream. 

rS 
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1 INTRODUCTION 

Molecular hydrogen (H2) is the most abundant molecule in the uni- 
verse and is closely linked to star formation via the star formation 
surface density rate - molecular gas surface density relation ( |Bigiel| 
|et al.|2008j l. Measuring rest-frame UV rotational transitions from 
the Lyman and Werner bands in absorption against a bright back- 
ground UV continuum is one of the few ways to directly measure 
H2 (see |Draine|201l[ for example). This technique probes diffuse 
gas with molecular fractions, /h2, of ^ to ~ 0.1 - denser 

molecular clouds are both dusty, and thus likely to extinguish UV 
light from a background source, and compact, such that there is 
a low probability of intersection with a sightline to a background 
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light source fHirashi ta et al.||2003[ |Zwaan & Prochaska|[2006| l. 
However, the lower molecular fraction systems that are detected 
give valuable insights into the environments and physical mecha- 
nisms necessary for the formation of H2. With this technique we 
can measure the physical properties of cool, dense gas over a large 
fraction of the age of the Universe, from the interstellar medium in 
the solar neighbourhood to proto-galaxies a few Gyr after the big 
bang. 

Since the initial detection towards the UV bright star ^ Per- 
sei jCamithers|1970^ , a large sample of sightlines exhibiting H2 in 
absorption from the Milky Way and its halo has been assembled. 
These observations have characterised H2 in diffuse molecular gas 
in the Milky Way plane ( [Savage et al.|1977) , the Magellanic clouds 
( [Tumlinson et al.|2002[|Welty et al.|2012^, high la titude sightlines 
out of the Milky Way plane jGillmon et al.|2006i |Wakker|2006l l, 
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in intermediate and high velocity clouds (IVCs and HVCs, Richter 



et al.|2003||1999|l, and in the Magellanic Stream jSembach et al 
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200 1[ [Richter et al.|2001[ l. A physical picture where H2 formation 
occurs pred ominantly on the surface of dust grains ( Shull & Beck-| 
with|l982 in clouds with total densities of n ~ 10 — 100 cm 
illuminated by the local UV radiation has been successful in re- 
producing both the observed H2 rotational population levels and 



molecular fractions in the Milky Way (for example Spitzer et al. 
|Jura|1975a| b) and the Magellanic clouds ( (Tumlinson et al. 
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20021. 



H2 has also been measured at redshifts 1.5 — 4.5, correspond- 
ing to lookback times of ~ 9—12 Gyr, in damped Lya (TVhi 
> 10^"-^ cm-^ DLA) and sub-damped Lya (10^'' cm"^ < A^hi < 
lO'^"'^ cm~^, sub-DLA) absorption systems seen towards bright 
background QSOs. In this redshift range, absorption features from 
H I and sometimes H2 are redshifted into the optical range, making 
them relatively easy to detect with large ground-based telescopes. 
The first unambiguous detection in a redshifted absorber was made 
by |Foltz et aLl ( [T988 |, and since then at least 16 further such sys- 
tems have been discovered (for example Ge & Bechtold| 1997 |^ Ge 
tt al. 2001 ; Levshakov et al. 2002 , |Cui et al.|2005[ |Ledoux e t al. 
2006„Noterdaeme et al.|2007| l. Approximately 10 per cent of DLAs 
have a molecular fraction /hj > 10~*'^, and these tend to be more 
metal rich and dustier JLedoux et al.'2003 i, and have higher veloc- 
ity widths ( |Noterdaeme et al. 2008 ) than DLAs without detectable 
H2. The inferred particle densities from the H2 rotational level pop- 
ulations and C I fine structure transitions are generally similar to 
those measured along local sightlines (~ 10 — 100 cm^''), but the 
ambient UV field, gas temperatures and gas pressures tend to be 
higher ^Hirashita & Ferrara 2005 , Srianand et al. 2005). 

No studies currently exist of H2 at lower redshifts, z < 1.5, 
outside the Milky Way halo. Until recently, the low number of DLA 
and sub-DLA systems known at low redshifts, together with the 
smaller light gathering power of spaced-based UV telescopes com- 
pared to large aperture ground-based optical telescopes, have made 
observing the Lyman- Werner bands in this redshift range imprac- 
tical. However, with the availability of the far-UV sensitive Cos- 
mic Origins Spectrograph (COS) on the Hubble Space Telescope 
(HST), molecular absorption can now be effectively detected for 
0.1<z< 0.8. 

In this paper we report the serendipitous detection of H2 in 
a sub-DLA at z — 0.56, the first such system analysed at a red- 
shift below 1.5 beyond the Milky Way halo. It has a high molecu- 
lar fraction given the total cloud neutral hydrogen column density, 
and we show that the associated metal absorption features seen re- 
quire the presence of three phases: a cold T ~ 100 K phase; a cool 
T ~ lO** K phase; and a warmer, probably coUisionally ionised 
phase. Based on the cloud properties we argue the absorber is likely 
caused by a tidally-stripped absorbing structure similar to the Mag- 
ellanic Stream embedded in a warm halo ~ 10 kpc from a nearby 
galaxy. 

The layout of the paper is as follows: Section 2 describes the 
data used; Section 3 describes how we identified lines and mea- 
sured the absorption line properties; and Section 4 describes the 
properties of the H2 absorption and the sub-DLA. We compare to 
theoretical models and discuss our results in Section 5, and sum- 
marise the main results of the paper in Section 6. When not ex- 
plicitly shown logarithms are to base 10, and we use a 7-year 
WMAP cosmology (Ho = 70.4 km s"^ Mpc-\ Q.m = 0.272, 
SIa = 0.728; |Komatsu et al.|20I T) where necessary. We use transi- 
tion wavelengths and oscillator strengths given by [Morton & Din-] 



Q 0107-0232 0l'^10™14.43= -02°16'57.6" 0.728 18.4 

Table 1. Properties of the background QSO towards which the sub-DLA 
is seen. Columns show the coordinates, emission redshift (measured from 
Mg II AA 2296, 2803 emission in the HIRES spectrum) and R band mag- 
nitude. 



Dataset 



Date obs. Exp. Time (s) Grating Ac (A) 



LB5H12010 


6 Nov 2010 


13905.728 


G160M 


1589 


LB5H13010 


18 Nov 2010 


13905.472 


G160M 


1589 


LB5H11010 


19 Nov 2010 


13905.376 


G160M 


1589 


LB5H14010 


24 Nov 2010 


13905.536 


G160M 


1623 


LB5H15010 


26 Nov 2010 


13905.504 


G160M 


1623 


LB5H16010 


7 Dec 2010 


13905.472 


G160M 


1623 



Table 2. Observations of Q 0107-0232 with HST/COS. Columns show the 
HST archive dataset name, the date observed, total exposure time, grating 
and central wavelength setting used. 



|erstein| p976) ; |Morton| ( (2003} , and [Vemer et al.] (1994}, and H2 
transition wavelengths from Bailly et al. ( [20I0[ >. 



2 DATA 

Transitions from the sub-DLA are measured in absorption against 
the continuum from the background QSO, Q 0107— 0232, at Zqso ~ 
0.728 (see Ta ble [T]l. This was di scovered by the Large Bright 
Quasar survey ( [Hewett et al.|I995] > and is one of a group of three 
bright QSOs with small angular separations on the sky. Spectra of 
these QSOs taken using the Faint Object Spectrograph (FOS) on 
the HST have been used to measure correlations in neutral hydro- 
gen absorption ( [Young et al.|2001 Petty et al.,2006) and in absorp- 
tion with galaxy positions jCrighton et al.|20I0 i across the three 
sightlines. 

Here we present higher resolution far UV spectra of 
Q 0107—0232 taken with the Cosmic Origins Spectrograph on 
the HST, and an optical spectrum taken with the High Resolution 
Echelle Spectrograph (HIRES) on Keck I. In our analysis we also 
make use of K band imaging of the QSO and archival UV FOS 
spectra. The FOS spectra were originally published bylYoung et al. 
(2001). We employ the combined spectrum used by[Crighton et al. 
(|2010jl, covering a wavelength range of 1572 — 2311 A at a typi- 
cal signal to noise (S/N) of 31 per 4 A resolution full width at half 
maximum intensity (FWHM). 



2.1 COS spectra reduction 

The COS spectra were obtained over a period from the 6th of 
November to the 7th of December 2010, as part of the Cycle 17 
proposal 11585. They represent a total exposure time of 23 hours 
across 30 orbits. Two central wavelength settings were taken with 
the G160M grating, each using 4 FP-POS positions to enable com- 
plete wavelength coverage from 1380 to 1850 A. Details of the 
exposures are given in Table|2] 

We used the CALCOS pipelin^to perform background sub- 



Version 2.13.6, [http : / /www, stsci ■ edu/hst/cos/[ 

[pipeline/[ 
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traction, wavelength calibration and extraction. The default back- 
ground extraction smoothing scale of 100 pixels resulted in poor 
background subtraction for our spectra, presumably because the 
pipeline was optimised for brighter targets. We found that chang- 
ing BWIDTH in the XTRACTAB calibration table from the default 
value of 100 to 20 significantly improved the background level such 
that the flux in strongly saturated features broader than the COS in- 
strument line spread profile was consistent with zero. 

Wavelength shifts are expected between visits and different 
wavelength settings due to temperature differences and uncertainty 
in the telescope pointing. The S/N in individual exposures is gen- 
erally too low (~ 2 per pixel) to reliably measure the centres of 
absorption features. Therefore we combined subsets of exposures 
grouping by FP-POS position, by visit (corresponding to a single 
dataset name in Table |2|, by grating central wavelength and by 
FUV segment to search for any shifts. Wavelength solutions were 
consistent across different visits and FP-POS values, but there are 
significant wavelength-dependent shifts between different central 
wavelength settings. To correct these, we measured the centroid for 
common narrow absorption features where two wavelength settings 
overlapped, and used these centres to calculate a wavelength offset 
as a function of position. We fitted these offsets with a linear depen- 
dence on wavelength, and then corrected for them such that FUV 
segment A Ac = 1627 matched the FUV segment A Ac = 1589 
setting, and FUV segment B Ac = 1589 matched the FUV segment 
B Ac = 1627 setting. The largest shifts applied in this way were 
0.1 A, corresponding to 20 km s^^, but they could result in a ~ 40 
km s^^ internal shift between the shortest and longest wavelengths 
of an exposure. These shifts are given in Table |Al] 

The scores of H2 absorption features distributed across the full 
spectral range enable a further check of the internal consistency 
of the wavelength solution. By measuring the centroid of these 
features and comparing to a single-component model of H2 ab- 
sorption, we discovered an additional wavelength-dependent shift 
(shown in Table |A2j. The magnitude of this shift is smaller 5 
km s~^) than that applied above, but still significant when fitting 
an absorption system with transitions spread across a large wave- 
length range. We removed this shift by subtracting a cubic spline 
fitted to the offsets as a function of wavelength position from the 
wavelength scale. 

To match the zero points of the COS and HIRES wavelength 
scales, we compared the N II 1084 and C I 945 features from the 
z = 0.56 system in the combined COS spectrum to their expected 
positions from the redshifts of the Fe II and Mg II lines from the 
same system measured in the HIRES spectrum. The wavelength 
zero point of the HIRES spectrum is known to better than 1 km s~^ 
relative to narrow Galactic Call absorption features seen in the 
spectrum. Both the N II and C I appear at redshifts expected from 
the HIRES Fe II and Mg II redshifts, and N II shows a similar com- 
ponent structure (albeit at the lower COS resolution). We conclude 
no correction to the wavelength zero point of the combined COS 
spectrum is necessary. 

We also measured the redshift of Galactic absorption features 
in the COS spectrum to confirm the zero point of the wavelength 
solution was correct. These are all saturated and possibly contain 
multiple components, so do not provide a stringent constraint on 
the zero point. However, they show no evidence of a systematic 
offset. 

After correcting each exposure for these wavelength shifts we 
made a combined spectrum in the following way. First we rebinned 
each exposure to a single wavelength scale with pixel width 0.0367 
A, ensuring Nyquist sampling. We used nearest-neighbour binning 



to preserve the spectra's noise properties, and checked that this did 
not introduce any significant wavelength shifts. The Icr uncertainty 
on each pixel was estimated empirically as the standard error on the 
mean of the contributing pixel fluxes. This is a slight overestimate 
of the true uncertainty, as the exposure times were not all identical. 
However, the uncertainties measured in this way are consistent with 
the standard deviation of the flux in regions free from absorption, 
and we believe this is a good estimate of the true uncertainty. 

Since the background level of the COS spectra is low, at small 
source count rates the flux distribution may be better described by 
Poisson rather than Gaussian statistics. However, in practice we 
find that for regions of our spectra with the lowest number of counts 
- the cores of saturated profiles - uncertainties in the background 
levels from the many contributing exposures makes a Gaussian flux 
distribution a good approximation. 

Finally we estimated the unabsorbed continuum level of the 
combined spectrum by fitting spline segments joining regions that 
appeared free from absorption. The resulting combined spectrum 
has a S/N of 10 per ^ 20 km s~^ resolution element at the contin- 
uum and covers a wavelength range from 1380 to 1850 A. 

2.2 HIRES spectra reduction 

The HIRES observations were performed on the night of 4th of Au- 
gust 2011. Four 1800 s exposures were taken using the red cross- 
disperser and a 0.861" width slit. Two wavelength settings were 
used to cover gaps in the detector. We used MAKEE to process 
each exposure, which subtracts the bias level and the sky back- 
ground, corrects for the echelle blaze, generates a wavelength solu- 
tion by identifying arc lines to yield a mapping from pixel number 
to wavelength for each echelle order, and extracts one-dimensional 
spectra for each echelle order. We then used custom-written Python 
code to coadd the individual orders for each exposure into a com- 
bined spectrum, and to infer the unabsorbed continuum level by 
fitting spline segments to regions free from absorption. The final 
combined spectrum has a S/N at 5000 A of 33 per 6.67 km s~^ res- 
olution FWHM, and covers a wavelength range 3890 to 8330 A. 

2.3 Imaging 

We acquired K band imaging of a 7' x 7' field around 
Q 0107-0232 using the High Acuity Wide field K-band Imager 
(HAWK-I) on the Very Large Telescope (VLT) during program 
383.A-0402. Five 180 s exposures were taken at four offset posi- 
tions on the 15th of September 2009. We used the HAWK-I pipeline 
recipes to process each exposure to remove the bias level and cor- 
rect for sensitivity variations using a flat-field. An astrometric solu- 
tion was measured for each exposure using SCAMP ( |Bertin|2006| l, 
then resampled to a common world coordinate system and coadded 
all the exposures with SWARP (Berlin et al.'2002). We determined 
the conversion between the measured counts and the magnitude by 
comparison to 2MASS magnitudes for objects in the field. The lim- 
iting magnitude reached is ~ 23.5 mag (AB) for a 3a detection of 
a point source. 



3 ANALYSIS 

3.1 Line identification 

Most of the transitions associated with the sub-DLA fall inside the 
Lya forest of the background QSO, and many are blended with ab- 
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sorption at different redshifts. We identified eacli absorption feature 
in the COS and FOS spectra in the following way. We first searched 
for Galactic absorption at the wavelengths of transitions typically 
seen in the Galactic interstellar medium (ISM; Sill A1526, CiV 
AA1548, 1550, Fell A1608, Ci A1657, and Mgi A2026/Znii 
A2026 were presenj^. Then we identified systems by the presence 
of either C iv (AA1548, 1550), O vi (AA1032, 1038), or H I Lya 
and Ly/3, starting at the emission redshift of the QSO and moving 
down in redshift to 2: = 0. Once these systems were identified, we 
searched for any further associated metal transitions such as Si IV, 
Si III, Si II, C III, C II. We found it was necessary to iterate this pro- 
cess several times, each time including line IDs from previous runs. 

The z = 0.56 sub-DLA was previously identified by Crighton] 
|et al.| ( |20I0^ by its many associated strong metal transitions in 
the FOS spectrum. Once we had made plausible identifications for 
lines at redshifts other than the z = 0.56 sub-DLA, we identified 
metals and molecular absorption lines from the Lyman and Werner 
bands for this system. Finally we assumed any remaining unidenti- 
fied absorption features were Lya. For this paper we focus on ab- 
sorption features associated with the z — 0.56 system. Absorbers 
at different redshifts are used only to identify blends with transi- 
tions from the sub-DLA. 



3.2 Kinematics and velocity structure of the sub-DLA 

H2 is expected to be found in gas with temperatures less than ~ 
5000 K - at higher temperatures molecules are destroyed through 
coUisional excitation ( |Shull & Beckw ith 1982). Therefore we ex- 
pect the H2 absorption features to be narrow, < 10 km s~^, and 
the COS spectra will not resolve the H2-bearing components. We 
use transitions covered by the higher resolution HIRES spectrum 
to inform us about the velocity structure of the absorbing gas, and 
then apply this to H2 and other transitions only present in the UV 
spectra. 

Figure^ shows the transitions at z — 0.56 detected in the 
HIRES spectrum: Mgll (AA2796, 2803), Mgl (A2853), Call 
(AA3934, 3969), and Fe II (AA2586, 2600). We also measure up- 
per limits on Al I, Fe I and Ca II. We fitted velocity components and 
column densities to these transitions using VPFIlr] The best-fitting 
values are given in Table [3] A single common velocity structure 
spanning ~ 200 km s~^ provides a good fit to all of these transi- 
tions, assuming line broadening is dominated by Gaussian turbulent 
motions rather than the gas temperature. The best fitting model is 
shown in Figure [T] Ca II and Mg I have the lowest ionisation en- 
ergies (11.87 and 7.65 eV respectively), and so a priori we might 
expect them to be associated with the cold environment where H2 
is found. However, the photoionisation analysis in Section [43] in- 
dicates that most of the Mg I and much of the Ca II probably arises 
in diffuse, photoionised gas distinct from the H2 . 

Component 6 has a Doppler width b (= V2a) of 20 km s^^, 
larger than is usually observed in low-ionisation metal transitions. 
This, together with the suggestion of correlated residuals in Mg II 
near the position of this component suggests it is in fact a blend of 
two or more narrower components. The quality of even the HIRES 
data is not sufficient to constrain the parameters of such heavily 

^ There is also absorption at the expected position of Zn I A2139, redwards 
of the QSO Lya emission. However, since this line has to our knowledge 
never been observed in the ISM of the Milky Way or redshifted DLAs, we 
leave the feature unidentified. 

http : / / www . ast . cam. ac . uk / ~r f c/vpf it . html 
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Table 3. Best-fitting Voigt profile parameters for each component in tran- 
sitions for the z = 0.56 system that are covered by the HIRES spectrum. 



blended components. However, as long as the distribution of un- 
resolved component widths is not strongly bimodal, the column 
density estimates for this component should be accurate jjenkinsi 
|I986[ l. We also measure A'^h2 independently of the velocity model 
assumed for H2 in Section [4/7] to ensure that the velocity model 
does not strongly bias our measurement of the molecular fraction. 

3.3 UV transitions for the sub-DLA 

We apply the Mgll velocity structure to models fitted to transi- 
tions observed in the lower resolution COS and FOS spectra. Using 
this velocity structure we were able to match the N II, Si II and O I 
profiles by varying the component line widths and column densi- 
ties. Several of the COS transitions that have measurable absorption 
and are not saturated or heavily blended with unrelated systems are 
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Figure 1. Observed transitions and upper limits at z = 0.56 in the HIRES spectrum. The thin smooth curve shows our best-fitting model using a single 
common velocity structure. The darker tick marks show components 5 & 6 that have associated H2. Points distributed around —0.25 show the residuals, 
defined as (flux — model) / (la error in flux). The two thin lines above and below this distribution mark the la deviation levels. Only upper limits are 
measured for Ca II and Fe I. 



shown in Figure [2] When fitting the COS spectra we use the tab- 
ulated line spread function provided by STScj^ linearly interpo- 
lated to the wavelength at the centre of each fitting region. We also 
measured column densities using the apparent optical depth (AOD) 
method (which assumes the transition is optically thin, [Savage &| 
|Sembach|1991| l, including a 5 per cent uncertainty in the contin- 
uum level. As the individual components are not resolved by the 
COS spectra, we quote these AOD measurements and give the total 
column densities for all components in aggregate. For transitions 
C I, N II, O I, O VI we were able to directly compare column densi- 
ties measured using both Voigt profile fitting and the AOD method. 
In each case they are consistent with one another. C II, C III and 
Si III are saturated, and lower limits are measured using the AOD 
method. The FOS spectrum provides an upper limit on A^siiv- Ta- 
ble|4]gives measurements and uncertainties, lower and upper limits 
for all of the transitions in the UV spectra. 

All H I Lyman series lines are saturated, and the Mg II velocity 
structure can reproduce the absorption well (but with larger, pre- 
sumably partly thermally broadened b values for each component). 
The damping wings measured at Lya in the FOS spectrum con- 
strain A*'hi = lO^^'^*'''^ cm~^, where the error is dominated by 
the systematic uncertainty in the continuum level (see Figure[3|. 



3.4 H2 velocity structure 

We measure H2 transitions from the J = — 3 rotational lev- 
els, with upper limits on J = 4 and 5. The asymmetric profiles 
for many of the H2 lines suggest there is more than one absorbing 
component. We were also unable to successfully fit the equivalent 
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Table 4. Total column densities for transitions in the z = 0.56 system 
observed in the COS and FOS spectra. A^hi is calculated from the damp- 
ing wings at hya in the FOS spectrum. C n, C III and Si III are saturated, 
and lower limits are calculated using the AOD method. The N I and Si IV 
values are 5a upper limits. The remaining values were calculated using the 
apparent optical depth of the transition with rest wavelength in the second 
column. Uncertainties given for these values are lu and include a 5 per cent 
uncertainty in the continuum level, which generally dominates the statistical 
uncertainty. ^: velocity models for Nil with a saturated central component 
allow higher column densities than this value and are still compatible with 
the data, but only by using a more complicated velocity structure than that 
fitted to the HIRES transitions. 



widths of the transitions using a curve of growth analysis with a sin- 
gle component. Therefore we fitted two H2 components, centred at 
the redshifts of the two central strong components at —26 km s~^ 
and km s^^ (components 5 and 6 in Tablejs} from the transitions 
observed in the HIRES spectrum. These are clearly separated in 
the resolution ~ 6 km s~^ HIRES spectra, but blended at the in- 
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Figure 2. Transitions for tlie z = 0.56 sub-DLA in the COS and FOS 
UV spectra. Tlie smootli red line shows Voigt profile models of the data. 
Absoiption that is not due to the named transition in each panel (usually 
from H2) is shown by dashed lines. Tick marks show the component posi- 
tions from Figure [T] dark ticks show components that have associated H2 
absorption. All transitions are covered by the COS spectra, apart from Si III, 
which is covered by the lower resolution FOS spectrum. 



strumental line profile of COS. However, the large number of tran- 
sitions over a range of oscillator strengths allow us to constrain 
velocity structure below the instrumental resolution. 

In Section [53] we show that the H2 is probably produced in 
a different environment to most of the metal lines. In this case, is 
it reasonable to use the component positions of the metals for the 
H2? It is likely that the H2 forms in the densest part of a cloud, 
and so will also be associated with the highest density (and thus 
strongest absorption) region of the warmer phase. So we believe 
this is indeed a reasonable assumption. 

We experimented with fitting the two components using VP- 
FIT, and found there were large degeneracies between the Doppler 
b parameter and the components' column density. One way to ro- 
bustly explore the b-N-z parameter space for the two components 
is to generate large grids of likelihood values as a function of the 
fitted parameters over plausible regions of parameter space. How- 
ever, in our case this proved to be prohibitively expensive compu- 
tationally. Instead we used a Monte-Carlo Markov Chain (MCMC) 
technique to sample parameter space. This samples parameter val- 
ues in proportion to the likelihood value at any point in parameter 
space. Thus from a set of initial parameter positions, a 'chain' of 
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Figure 3. Constraints on Afjji from the Lya transition. The histogram 
shows the data, the dashed line the continuum, and the lower green points 
the residuals as defined in the caption of Figure [T] The thick red solid 
curve shows the best-fitting A^hi = 10^^'^ cm"'^ for a single component 
with b = 20 km s~^, and the thinner upper and lower solid curves show 
Nm = lO"-'* and lO^^'^cm-2. 



Component 5 



Component 6 



J 


logioA^(cm-2) 


b(kms-l) 


logiQ Af (cm^2) 6(kins~i) 





16.17 ± 0.25 


7.1 ± 1.1 


15.63 ±0.39 4.3 ±0.7 


1 


17.05 ±0.28 




16.42 ±0.40 


2 


16.19 ±0.19 




15.65 ±0.25 


3 


15.77 ± 0.12 




15.47 ± 0.18 


4 


< 14.5 




< 14.5 


5 


< 14.5 




< 14.3 



Table 5. Column densities for the two components showing H2 based on 
Monte Carlo Markov Chain fitting. The en'ors are Icr. The J = 4 & 5 
values are 5cr upper limits. 



parameter values is generated by a stochastic walk through param- 
eter space with distributions approximating the Bayesian posterior 
probability for each parameter. 

We generated posterior parameter distributions using the 
package EMCEE (The MCMC Hammer, ,Foreman-Mackey et al.| 
|2012[ l. We fitted for each component's redshift and b parameter, 
and for the column density of each rotational level using the 56 
transitions shown in Figure]?] All transitions for a component were 
constrained to have the same b value, which is generally observed 
to be the case in local sites of H2 absorption for at least J levels 
< 3 ( |Spitzer et al.|1974) . Fitting transitions from each rotational 
level individually also gives b parameters consistent with a single 
value. Even after correcting the wavelength scale, residual wave- 
length shifts of ~ 1 km s^^ remain, so we also allowed a small 
wavelength shift for each fitted region. 

Table]5]gives the parameter estimates and la errors for the H2 
column densities and b parameters for each rotational level. The Icr 
regions are determined by marginalising over all other parameters 
and finding the narrowest region that encompasses 68.3 per cent of 
the samples. We choose the parameter values to be at the centre of 
these la regions. The absorption model with the set of parameters 
that maximises the likelihood is shown in Figure]?] 
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Figure 4. The H2 absorption model generated using parameters that maximise the likelihood of the data. The histogram shows the flux divided by the 
continuum (the darker portions were used to calculate the likelihood), the smooth curves show the model. The residuals - (data - model) I la - ai'e shown 
centred on —0.25, scaled such that the light gray lines above and below —0.25 are the ztlcr range. A two component model with a single b-pai'ameter for each 
component and redshifts corresponding to components 5 & 6 in Table|3]reproduces the data well. 



4 ABSORPTION SYSTEM PROPERTIES 
4.1 Metallicity 

The metallicity, Z, can be estimated from an element X using the 
log of the ratio of the abundance of element X in the absorber, 
Nx/Nh, to the solar abundance 

(iVx/iVH)obB 



[X/H] = logi 



(1) 



'1° {Nx/Nh)q • 
Due to a charge transfer between O and H we expect the ratio of 



the number densities no: /no and nm/nu to be the same (Field 
|& Steigman||197l) , provided the majority of O is in the form of 
O I and 11. In the presence of many high energy ionising photons, 
this is no longer true due to different absorbing cross sections of O I 
and H I (see |Prochter et al.|2010j l, but the absence of Si IV argues 
against a hard radiation field for this system, and the best fitting 
Cloudy models do not predict significant amounts of O in higher 
ionisation states (O VI is seen, but we argue this occurs in a hotter, 
coUisionally-ionised phase). Oxygen also shows little depletion (< 



0.3 dex) onto dust grains across a range of environments in the 
ISM of the Milky Way ( ,Jenkins|2009^ , so should provide a good 
estimate of the metallicity. 

We find [O I/H I] = -0.86±0.24, or ~ 0.14 Zq. For the pho- 
toionisation analysis in Section [43] we assume the metallicity is the 
same across the entire complex. We caution that in one of the few 
cases where the metallicity has been measured for individual com- 
ponents in a single absorption system, metallicity differences of a 
factor of ten have been observed l |Prochter et al.|20To !l. However, 
given that the dispersion in A^Mgii/A^Feii across the system is not 
excessively large (the largest log difference between components is 
0.38), the assumption of a constant metallicity seems reasonable. In 
particular, the two components with H2 do not show significantly 
different ion abundance ratios compared to the entire system. 



4.2 Dust 

|Noterdaeme et al.| ( |2008^ have found a correlation between the pres- 
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ence of dust and the likelihood of observing H2 in DLAs, con- 
sistent with the main formation mechanism for H2 being on the 
surface of dust grains. We can measure the dust content by com- 
paring elements known to deplete strongly onto dust grains (Fe, 
Mg) to those with low depletion (O). We assume negligible ion- 
isation corrections, but applying corrections from the best-fitting 
Cloudy model in Section [43] does not change our conclusions. 
We find [Fe/0] = -0.15 ± 0.08 and [Mg/O] = -0.25 ± 0.07 
for the entire system, indicating mild dust depletion. If we also as- 
sume solar abundance ratios we can estimate the dust-to-gas ratio 
normalised by the value in the solar neighbourhood as 

K = 10[^/H1(1-10P^/^1), (2) 

where X is an element that does not deplete strongly onto dust (for a 
derivation of this expression see the appendix of |Wolfe et al.|2003] l. 
Using oxygen gives logjg k — — 1.4lg g, which is at the lower 
end of the range of values found in other higher redshift systems 
showing H2 jLedoux et al.|2003) . 

4.3 Temperature constraints from line widths 

The linewidths of absorption components in the HIRES spectra 
can be used to constrain the temperature of the gas using the re- 
lation b = ^2kT/m, where m is the mass of the ion and T is 
the temperature. This constraint is an upper limit, as there can be 
large-scale turbulent motions in addition to thermal broadening, or 
the line may not be resolved. Indeed, we fitted each component 
in Mgll, Fell, and Call with a single b parameter value across 
all three transitions, consistent with turbulent broadening dominat- 
ing over thermal broadening. Due to the relatively large masses for 
these elements, only component 2 gives a constraining upper limit 
of 6 000 K, typical of temperatures in the warm neutral medium 
in the Milky Way ISM. The H2 linewidths give upper limits to the 
temperature of 8 800 K and 5 200 K for components 5 and 6, but we 
argue below that the physical conditions in the H2 gas are probably 
different to those of the gas where most of the metal lines arise. It 
is also possible that the H2 widths are substantially broadened due 
to turbulent motions. 

In conclusion, there are no strong temperature constraints 
from the linewidths. The relative column densities of the H2 ro- 
tational levels provide an independent measure of the temperature, 
discussed in the next section. 



4.4 H2 excitation temperature 

The ratios of H2 column densities in different rotational levels can 
be expressed as excitation temperatures, assuming a Boltzmann 
distribution across the levels (see |Draine(201 l[ l: 

AT ^ ™P ^ 



J 



1 



(3) 



Here Nj is the column density for molecules in rotational state J, 
and gj = (2 J + 1)(2/ + 1), where / = if J is odd or 1 if J is 
even, is the statistical weight of J. = 85.36 K and Tjo is the 
excitation temperature from level J to J = 0. 

Fig.|5]shows an excitation diagram for the column densities of 
the J = — 3 transitions for the two H2 components. For densities 
above nu ~ 100 cm^'^ the J = and J — 1 levels are domi- 
nated by collisional processes ( [Dalgarno et al.|1973^ , and thus Tio 
represents the kinetic temperature of the gas. The uncertainties on 
the column densities for the J = and J = 1 levels yield a lower 
limit on Tio for each component at la (2a) limits of 123 K (64 K) 



• Component 5 

T = 2iOK 

▲ Component 6 
- T = 290A' 
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Figure 5. Excitation diagram for the two H2 components. The shaded re- 
gions show the 1(T error range. The slope of a line joining each pair of 
points is inversely proportional to the excitation temperature between those 
two J levels. An illustrative temperature is shown for each component, but 
the population levels are not expected to follow a Boltzmann distribution 
characterised by a single temperature. 



for component 5 and 77 K (37 K) for component 6. We also show 
two illustrative temperatures corresponding to the populations for 
J = — 3 for each component. However, we caution that differ- 
ent physical processes affect the populations of these levels |Jura| 
|1975b| l, so it is not expected that a single temperature should match 
all four levels. 



4.5 Cloudy modelling 

In this section we attempt to generate a simple single-cloud model 
illuminated by a UV radiation field that can reproduce all the ob- 
served column densities. We compare to the total column densities 
for all components, since the individual component columns are 
not well constrained for the O I, C I, and Si II transitions or the sat- 
urated transitions (C II, C III and Si III). Given the large range of 
transitions present with widely differing ionisation energies, it is 
likely that there are several different phases present, and a single 
cloud model is unlikely to be able to reproduce all the observed 
species. Below we find that a single model can reproduce the ma- 
jority of the low-ionisation metal transitions, but Section [53] shows 
that multiple phases with different densities and temperatures are 
required to explain all the absorption. 

We use models generated with version 8.01 of CLOUDY, last 
described by [Ferland et al.| jl998[ >, to estimate the physical con- 
ditions in the absorption system. All models assume solar abun- 
dance ratios, constant gas density, and an absorbing geometry of a 
thin slab illuminated on one side by an incident radiation field per- 
pendicular to the slab surface. The radiation field includes the cos- 
mic microwave background at the redshift of the absorber. We then 
compare four scenarios: a cloud in an intergalactic medium-like 
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(IGM) environment, in an ISM-Iike environment, close to a star- 
burst galaxy, and illuminated by an AGN-dominated spectrum. The 
IGM-like model is free of dust with a radiation field given by the 
UV background spectrum from |Haardt & Madau| ( |2012l l, including 
contributions from quasars and star-forming galaxies at the redshift 
of the absorber. It has a radiation field strength at 912 A, J^jxiM ~ 
6.08 X 10"^^erg s"^cm"2Hz"^sr"\ The ISM-like models have a 
radiation field similar to the Galactic ISM (J^^'^ ~ 400JJ^\gm). 
which is dominated at UV wavelengths by the spectral shape of 
hot stars, and a dust grain composition similar to that measured in 
the ISM. Even though the ISM-like models use solar relative abun- 
dances, the gas phase abundance ratios are substantially different 
from solar due to differential depletion of metals onto grains. The 
starburst CLOUDY models assume the absorbing cloud is lOkpc 
from the galaxy and an escape fraction for UV light of 3 per cent, 
in addition to the IGM-like radiation field described above, without 
dust grains. They have jj*^^ ~ 2700Jj^\GM- 

The starburst galaxy 

spect rum used was generated using STARBURST99 l |Leitherer et al.| 
1999[ l for a star formation rate of 20 Mq yr"^ The AGN models 
use the default tabulated AGN spectrum from CLOUDY with a nor- 
malisation Jj'^^ ~ 3000Jj^\?3M 

and do not include dust grains. 

For each scenario we generate a grid of models for a range of 
ionisation parameters, metallicities and total A^hi- We estimate the 
ionisation parameter U, defined as the ratio of the densities of ion- 
ising photons to hydrogen atoms, using the observed total column 
density ratios iVMgi/A''Mgii, Nm/Nmi, and A''siii/A'"siiii- Using 
ratios of ionisation states for the same element avoids any effects 
that might alter the column densities of ions for different elements 
in different ways, such as non-solar abundance ratios or differen- 
tial dust depletion. We generate the likelihood of each parameter 
(U, Z and A''hi) for the grid of models based on the observed ra- 
tios, and include a Gaussian prior on the metallicity centred on the 
[O I/H I] metallicity measurement with width a equal to the la un- 
certainty on the metallicity. For all three scenarios, only a relatively 
narrow range of U values correctly reproduces the observed ratios. 
The likelihoods are only weakly dependent on the total A'^hi; we 
assume A'hi = 10^^'^ cm""^, which results in models that best re- 
produce the observed metal column densities. 

Once we have found the most likely U value, we compare the 
predicted column densities to the observed transitions with mea- 
surements or limits, and assess which scenario reproduces the ob- 
servations best. We first compare to the AGN models. These are 
the only models with a hard enough spectrum to produce sufficient 
A^ovi to match the observed value. However, at the same time they 
overpredict the amount of Si IV, Fell, Mgl and Mgll by one or 
more orders of magnitude. Thus it is more likely the O VI arises in 
a coUisionally-ionised phase separate from the low-ionisation tran- 
sitions, as is observed in other systems (for example Riba udo et al.| 
|201 Ij l, and we do not compare to the high-ionisation species (Si IV, 
O VI) for the remaining models. 

The most likely model for the ISM scenario underpredicts the 
abundances of Fe II by more than an order of magnitude, and Ca II 
by many orders of magnitude, as both are expected to be heav- 
ily depleted onto dust grains. This suggests the depletion pattern 
in the z = 0.56 sub-DLA is different from that in the Milky Way 
ISM. This model also underpredicts C I. The starburst scenario fails 
to reproduce the observed Nt^i/Nt^u, and also severely underpre- 
dicts Ca II, Fe II and C I. The IGM-like model gives the best fit to 
the observed data, and its predictions along with observed column 
densities are shown in Figures [6] & [7] Figure |6] shows that col- 
umn densities for O I, Mg I, Mg II, Fe II, N I and N II are reasonably 
well matched. The remaining small deviations from the predictions 
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Figure 6. Cloudy predictions (lines) and observed column densities 
(points) for Mgl, MgO, Nl, Nil, Ol and Fell for the IGM-like model 
with Z = 0.13 Zq. Points are plotted at the most likely U values based 
on constraints from the A^Mgl/^Mgll, ^Nl/N-f^n, and A'^siii/Afgim ra- 
tios (Fe II is offset for clarity). Filled points coiTespond to solid lines of the 
same colour, open points to dashed lines. Errors on the column densities are 
smaller than the marker sizes. Given the uncertainties in the shape of the 
ionising spectrum and uncertain relative abundances, the predicted column 
densities are veiy close to the observed values. 



could be due to a slightly different incident UV continuum from the 
one assumed, or enhanced or depleted elemental abundances rela- 
tive to the solar values assumed. For example, 0.2 dex less A'nii is 
observed than is predicted. This may be due to a nitrogen under- 
abundance, often observed in similar A^hi systems at low redshifts 
(for example Battisti et al. 20121 and in DLAs at higher redshifts 
( [Pettini et al. 2002, Prochaska et al.,2002^ . Figure|7]shows that Si II, 
Si III, C II and C III columns are reproduced well. However, there 
is still 0.5 dex too little A'^caii and 1 dex too little A'^ci predicted. 
In Section [53] we suggest a scenario to explain this discrepancy 
between the models and observations. 

We also ran CLOUDY models using constant pressure clouds 
instead of constant density. The motivation for these was to simulta- 
neously include cool, lower density gas at the edge of the cloud, and 
higher density, cold ^ 100 K gas at the core of the cloud where H2 
can survive. In these models we also included contributions from 
cosmic rays, which can be important for cold molecular regions. 
Although significant amounts of H2 can co-exist with many of the 
metal transitions observed for these models, they still cannot cor- 
rectly reproduce the Ca II or C I columns. 



4.6 C II fine structure absorption 

Singly ionised carbon (C^) has electronic structure 2s^2s^2p 
where the outer shell has a configuration ^Pj, and thus fine struc- 
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Figure 7. Cloudy predictions (lines) and observed column densities 
(points) for Cl, C II, C III, Si II, Si III and Ca II for the IGM-like model with 
Z = 0.13 Zq. Points are plotted at the most likely V values based on con- 
straints from the Nyigl/l^Mgll, ^TSu/^Nll, and A^sillZ-'Vsilll ratios (Si 
and C III points are offset for clarity). Filled points coirespond to solid lines 
of the same colour, open points to dashed lines, and open points with dotted 
edges to dotted lines. Where error bars are not shown they are smaller than 
the marker size. With the exception of Cl and Call, the column densities 
are reproduced well. As discussed in Section [53] the excess Nq j is likely 
due to C I being in a self-shielded cold phase where the H2 and possibly 
part of the Ca II resides. 



ture splitting occurs between the J = 1/2 and J = 3/2 levels. The 
transition from the ground state to these levels produces C II and 
ClI* absorption respectively. 

The ratio of ClI* to ClI column densities has been used 
to estimate the star formation rate inferred from the cooling rate 
in DLAs at high redshift ( | Wolfe et all] [2003} . We find an up- 
per limit on ClI* from A1036 of A^cii* < lO"'^ cm'^ and 
A^cii > 10^'' '^ cm" ^. Assuming constant n(C II* )/n(C II) over 
the entire complex, we find the ratio A^'cii* /Ncii < 0.8, consis- 
tent with ratios measured in higher redshift DLAs (for example!Sri-| 
[anand et al.| 2005[ l and local environments. Following the assump- 
tions described by |Morris et al.i (1986), we can estimate the electron 



density in cm 

iVcii 



' using the expression: 



iVcii 



3.9 X 10" [1 + (0.22- 



(4) 



We use Up/ne corresponding to the ionised H fraction of 0.70 from 
the best-fitting CLOUDY ionisation model to find rie < 10 cm"''. 
Using A'hi we can estimate the thickness of the absorbing cloud 
as Nm/rie - This gives a lower limit on the cloud size of 1 pc. This 
limit is not necessarily related to the density or size of the H2 gas, 
as we argue in the discussion that most of the C II is due to a warm 
ionised phase separate from the H2 . 



4.7 Molecular fraction 

The molecular mass fraction is estimated by 

/h, =2iVH,/(7VHi + 2iVH,), 



(5) 



assuming most of the hydrogen associated with the H2 is neu- 
tral. In this case, as for many other QSO absorption systems, it 
is not clear how to divide the total TVhi measured from the damp- 
ing wings between different absorbing components, and in prin- 
ciple each A^H2 component could have a different /hj value. To 
calculate /hj we use the total from both components and 
conservatively assume all A^'hi is associated with the H2, mean- 
ing is effectively a lower limit. This gives a molecular fraction 
of logjg ~ —1.93 ± 0.36. As we discuss in the next section, 
given the total A^hi, this is an usually high molecular fraction com- 
pared to most other higher redshift systems and sightlines in the 
Local Group. Therefore we may be concerned that a different ve- 
locity model to the one we have used permits a much lower /hj . To 
calculate a lower fimit on independent of the velocity model, 
we measure the column density of the lowest oscillator strength 
transition available for each rotational level (J — 0, A1108.1; 
J = 1, A1008.5; J = 2, A934.1 and J = 3, A952.3) using the 
AOD method. This gives a lower limit of Nii2= 10^^'^ cm"^ or 
/h2 > 10"^, again assuming all of the A^'hi is associated with the 
Nu2- This is still a high value relative to local H2 systems with 
similar A^'hi. 



5 DISCUSSION 

5.1 Physical conditions in the H2 cloud 

To consider this system in the context of other H2 detections in ab- 
sorption, we plot for local and higher redshift H2 sightlines as 
a function of the total hydrogen column density and in Figure 
[8] It is apparent that the z — 0.56 sub-DLA (solid circle) has an 
unusually high /hj given its A'^hi compared to sightlines through 
the plane of the Milky Way (cyan inverted triangles), or through the 
Magellanic clouds (green squares and red diamonds). 

Before we discuss the likely origin of the z = 0.56 sub-DLA, 
we examine the physics underlying the /h2 distribution as a func- 
tion of Nil and 7Vh2 ■ The left panel shows a clear bimodality in 
the /h2 —Nh distribution between high A'h, high /h2 sightlines 
at the top right, and lower /h2 sightlines, generally with much 
lower A'^H. The right panel shows this is actually a bimodality be- 
tween A'h2 values < 10^® cm"'^ and > 10^* cm"^. This can be 
understood as the onset of H2 self-shielding against UV dissociat- 
ing photons (for example [Hirashita & Ferrara|2005 Gillmo n et al.| 
|2006[ l. An analytic approximation from Draine & Bertoldi ( 1996]l 
shows that > 97 per cent of H2 -dissociating photons are blocked 
by self-shielding once A''h2 ~ 10^^ cm"^. Once H2 becomes self- 
shieldetj^ the dissociation rate drops and H2 accumulates rapidly to 
the formation-destruction equilibrium value predicted by the mod- 
els by |McKee & Krumholz|p010^ . These models are shown at the 
top right in each panel for two metallicities; solar and Z — 0.2Zq, 
the metallicity of the SMC. They were calculated using equations 
4, 5, 7 and 8 from |Kuhlen et al.] ( |2012i and assume the ISM is in a 
two phase equilibrium between a cold neutral medium and a warm 
neutral medium (for example [Wolfire et al.|1995) . The solar metal- 
licity model reproduces the mean /h2 for the Milky Way sightlines 



^ Dust shielding only becomes important at total H columns of 
10'^^ cm"^, assuming solar metallicity. 




Figure 8. The molecular fraction versus the total A'^jj . Compared to measurements in the Magellanic clouds jTumlinson et al.|2002 



jTum li 

the disk of the Milky Way ^Savage et al.|1977) , in intermediate- (Richter et al. 2003) and high-velocity clouds (Richter et al. 1999 , and in most other QSO 
absorption systems at z > 1.5 (from the compilation by Noterdaeme et al.|2008| With additions from |Petitjean et al.,2002nReimeirs"et al.|2003]|Tumlinson| 
|et al.|20i0| and [jorgenson et al.|201Q ) the z = 0.56 H2 system has an unusually large /hj for its total H column. The local systems that appear most similar 
to this absorber are seen along sightlines through the Magellanic Stream ^Sembach et al.|2001 [[Richter et al.|200H . Error bars at the comer of each plot show 
the typical uncertainties on fn^ , A^tot and A'^H2 ■ Shading shows the regions in which H2 cannot be detected for a limit of Nn^ < lO^*'^ cm~^, a typical 
threshold for the spectra used for the different surveys. The thin lines at the top right in each panel show analytic models from Krumholz and McKee for two 
illustrative metallicities. The three thin lines at the lower left are the analytic models described in Section [5T] 



[Welty etal.[20T2t , along 



through shielded H2 regions reasonably well, although |Welty et aL] 
( [2012 ') point out these models overpredict the A^h at which this 
transition occurs. 

The A'^H at which there is sufficient shielding from the UV 
field to form large amounts of H2 varies depending on the dust to 
gas ratio, the strength of the UV field, and the H2 linewidth, so the 
A'h at which the transition from optically thin to optically thick 
occurs can change froin sightline to sightline. In the plane of the 
MW disk, the transition from low to higher /hj takes place around 
A^H = lO^"-^ - 10^^ cm"^ It occurs at higher Nh in the LMC and 
SMC, both because their lower metallicities (0.5 and 0.2Zq re- 
spectively, see the appendices from |Welty et al.|1997[|199"9| l result 
in a lower H2 formation rate on grains, and due to an increased UV 
field compared to the Milky Way ( [Tumlinson et al.|2002[ l. [Grilmon| 
|et al.| l |2006l l have also shown the large variation in along differ- 
ent sightlines implies that each sightline intersects a small number 
of molecule-bearing clouds. 

In addition to comparing to the two-phase equilibrium models 
of McKee & Krumholz, we can compare to simple analytic models 
that apply to diffuse H2 in the partially shielded regimes. Following 
the Appendix from jjorgenson et al.j j2010^ , the dissociation rate 
in s~^ due to photons with energies corresponding to the Lyman- 
Werner bands is given by 



1.1 X W^TTjlfSMM, 



(6) 



shielding, which can be calculated using the analytic expressions 
from Draine & Bertoldi| ( fT996) and |Hirashita & Ferrara|p005t . In 
formation-dissociation equilibrium, the molecular fraction can be 
approximated by 



/h2 =2 



(7) 



where TZ is the formation rate of H2 on dust grains in cm^ s~^, 
riHi is the Hi particle density in cm~^ and k is the dust to gas 
ratio relative to that in the solar neighbourhood as defined in Sec- 
tion |4.2| Note that the formation rate term used by jHirashita & Fer^ 
[rara (2005 1, -Rdust, is equal to nTZ. Rearranging these expressions, 
we can estimate the particle density in the cloud as: 



nm = 74 cm 



/h2 
0.01 



•J u 



tLW 

•^v, SN 



5si 



0.01 



(8) 



where j]^^-^ = 3.2 x 10~^° erg s"^ cm"^ Hz"^ sr"^ and 7^sN = 
3 X 10~^^ cm'^ s~^ are typical values measured in the solar neigh- 
bourhood ( |Habing| 1 968"| j Jura| 1 974^ . 

The three curves at the lower left of each panel in Figure [8] 
show the molecular fractions estimated with equation (7) for three 
illustrative combinations of nTZnm/ Riisi- The upper curve and 



where jj;* is the strength of the incident radiation field in 
erg s'^cm^^Hz^^sr^^ at 1000 A, and 1 — SshicU is the fraction 
of Lyman- Werner photons processed by dust or scattered due to H2 



^'^ and K = 1, 0.337?.sN, 2JJ;^n respectively. The lower 



10 Jj,, SN 

curve has k — 0.04, nni ~ 5 cm"'', 7?.sn and 0.01J^*sn- There 
fore the observed variation in can be explained by reasonable 
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variations in the combination of UV field strength, particle den- 
sity and H2 dust formation rate. If the z — 0.56 system is in H2 
formation-dissociation equilibrium, the combination of low A^hi 
and high molecular fraction suggests that it is either in a weaker UV 
field, has an increased H2 formation rate, a higher nn compared to 
the solar neighbourhood, or some combination of these three. We 
can use limits on the column den sities of the J — 4 and 5 levels 
to put upper limits on jj;* ( Jura"'l975b). These upper limits are 
not very stringent, however, due to weak limits on the column den- 
sities, and constrain J^f < 10^^** erg s^^cm~^Hz^^sr^^. This 
is consistent with the different UV background values assumed in 

10"^^ for the IGM-like 
for the ISM-like see- 



the Cloudy models, which range from ' 
scenario to ^ 10^'^" erg s^^cm^'^Hz^^sr 
nario. 

Using and measured in the z — 0.56 sub-DLA, 
the measured metallicity and equation (8), we find densities of ~ 
1 — 4 cm''' for a UV background incident radiation field, and ~ 
70-480 cm~^ for a Milky Way ISM-like radiation field. The lower 
density range corresponds to cloud thicknesses of ~ 3 — 10 pc, the 
high density range to ~ 0.002 — 0.15 pc. The upper end of this 
density range is consistent with values measured for higher redshift 
H2 systems using C I fine structure transitions, but would result in 
extremely small cloud sizes. 

The low total column density of the z — 0.56 system sug- 
gests that it does not pass through the ISM of a galaxy. Returning 
to Figure. [5] we see that local systems with similarly low total A'^h 
and almost as high are sightlines through a high velocity cloud 
([Richter et al.||I999|l and the Magellanic Stream ( [Sembach et al.| 
|2001[|Richter et al.|2001^ . These clouds have sub-solar metallici- 
ties (0.3 — 0.5 solar), and are most likely tidally stripped from the 
Magellanic Clouds (for the Magellanic Stream) or the Milky Way. 
[Sembach et al.| ( [2b01| ) estimate the density of the H2 -bearing cloud 
they observed in the Magellanic Stream to be 0.3 — 3 cm~'^ with a 
photoionisation rate at least a factor of 10 smaller than the Milky 
Way ISM value. The H2 formation timescale for these low densi- 
ties is around 1 Gyr, a large fraction of the estimated lifetime of 
the Magellanic Stream (for example |Besla et al.|20I0[ >. Therefore 
they favour a scenario where H2 is not formed in place, but has 
survived the tidal stripping process and persists due to a combina- 
tion of self-shielding and the lower ambient UV field compared to 
the LMC ISM. Such a scenario could also be responsible for the 
z — 0.56 absorber. 



5.2 Connection to galaxies 

The K band imaging around Q 0107—0232 has a seeing FWHM of 
0.8", and shows two possible galaxy candidates less than 1.2" from 
the QSO sightline. Figure [9] shows a 5"x5" region centred on the 
QSO. The QSO image has been subtracted using the point spread 
function of a nearby star. The two galaxy candidates are seen to 
the North-West (Gl) and South-West (G2). Assuming they are at 
the redshift of the absorber, they have luminosities of 0.7L* (Gl) 
and 2L* (G2), and impact parameters of 10 kpc (Gl) and 11 kpc 
(G2), both smaller than the median impact parameter of 33 kpc 
for galaxies associated with sub-DLAs found by |Rao et al.|p01l] ). 
Therefore it is likely that at least one is associated with the absorber, 
on scales typical of the separations between the Milky Way and 
high-velocity clouds (10-60 kpc, see |Putman et al.||2012 | and the 
references therein). Clearly it is desirable to measure redshifts for 
these candidates. 
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Figure 9. A K-band image of the QSO showing two nearby galaxy can- 
didates, Gl and G2. The QSO image has been subtracted using the point 
spread function of a nearby star. The region near the centre of the QSO still 
shows significant residuals and has been greyed out. If Gl and G2 are at the 
redshift of the QSO then they have impact parameters of 10 kpc and 11 kpc, 
and luminosities of 0.7 L* and 2L* respectively. 



5.3 Three different gas phases in the sub-DLA 

Figures [6] &|7] show the total hydrogen particle density for the ma- 
jority of metals observed in this system corresponding to the ioni- 
sation parameter assuming the normalisations of the incident radi- 
ation fields are correct. The most likely model corresponds to hy- 
drogen densities from 10~^ to 10~^ cm~^. Even assuming a factor 
of ten uncertainty in the radiation field strength, this is much lower 
than the typical densities where H2 is seen in both our galaxy and 
in other H2-bearing DLAs (nn = 10 — 100 cm'^). Therefore the 
gas traced by most of the metal absorption is probably in a different 
environment to that in which the H2 resides. This is also likely the 
cause of the excess A^'ci over that predicted by the CLOUDY mod- 
els. C I is often seen in dense components showing H2 (for exam- 
ple Srianand et al. 2005 1 and can have extremely narrow linewidths 
corresponding to temperatures of ~ 100 K ( [Jorgenson et al.|2009[ 
[Carswell et aI.|20rT| , indicating it occurs in the same environment 
as H2. Thus most of the Cl and some Call may be from a high 
density region co-spatial with the H2 . 

As discussed in Section [43| the presence of O VI is unlikely to 
be explained by photoionisation by a hard UV field, and indicates 
the presence of a hotter medium. We conclude that the absorption 
is due to gas in three phases: a photoionised cool phase at ~ 10* K 
in which most of the metal transitions we see are produced, a cold 
^ 100 K phase where the H2 and C I absorption occurs, and a hot- 
ter phase where O VI is produced. The H I column is likely split 
between the two cooler phases. A similar multi-phase environment 
is also seen in other hi gher redshift sub-DLAs that show molecular 
absorption ( Milutinovic et al. 2010) ^] 

The Magellanic Stream and most other HVCs comprise 10* K 



Milutinovic et al. did not report an Ovi detection, but as this system 
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ionised gas that is seen in Ha emission, T > 10° K hot gas pro- 
ducing Ovi absorption, and cold neutral gas with H2 ( [Sembachl 
|et al.|20()3{|Fox et aLpOlO ). Taken together, the existence of these 
three phases, the high molecular fraction with a low total column 
density, and the proximity of a possible ~ L* galaxy suggest the 
z = 0.56 absorber is due to a tidally stripped feature analogous to 
the Magellanic Stream. 

5.4 Incidence rate of H2 in low redshift sub-DLAs 

Due to the need for bright targets observable with space-based UV 
spectroscopy and their low incidence rate, very few DLAs and sub- 
DLAs have been found at low redshift. Until recently only ~ 10 
DLAs at redshifts < 1 were known, and only a handful of these 
have coverage of H2 Lyman-Werner bands. With the availability of 
COS, the number of such systems is being increased dramatically, 
and due to its far UV wavelength coverage the presence of H2 can 
be easily detected. 

[Battisti et al.| ( |20i2) l present a sample of 2 DLAs and 6 sub- 
DLAs at z < 0.35, serendipitously discovered along sightlines as 
part of a large COS program. Like the sub-DLA presented here, 
they were not pre-selected by the strength of their metal lines or 
other properties that might influence the likelihood of detecting 
molecules. Interestingly, they also discovered a sub-damped sys- 
tem with H2 absorption at 2; = 0.2477. Taking this sample to- 
gether with the system in this paper and assuming binomial statis- 
tics, we find the expected incidence rate of DLAs and sub-DLAs 
showing molecular hydrogen at 10^* cm~^ at low redshift 

to be 2/9 = 22 per cent (with a 95 per cent confidence level lower 
limit of 4 per cent), rising to 33 (5) per cent if we consider only the 
sub-damped systems with A^hi < 10^" cm~^. This is a surprisingly 
large fraction given that sub-DLAs are often found to be highly 
ionised absorbers with < 10 per cent of their hydrogen in the form 
of Hi. 

If we think that the absorption cross section for H2 is dom- 
inated by cold gas associated with Local Group-type systems (the 
Magellanic Stream for example), then this would be consistent with 
this high incidence rate. Richter (2012) shows that one can explain 
30-100 per cent of the observed incidence rate of systems with 
A^'hi > lO^'^'^ cm~^ as intermediate- and high-velocity clouds dis- 
tributed around galaxies with Hi masses between 10*'^ and 10^" 
Mq in a similar way as is seen around M3 1 and the Milky Way. 
As discussed in the previous section, some HVCs also show rela- 
tively high molecular fractions, and in terms of A^'hi and /h2 HVCs 
are the local systems most analogous to the system analysed in this 
paper. 

It would be interesting to perform a systematic search for H2 
in further 10^" cm"^ < A'^hi < lO^'^'^ cm~^ sub-DLAs at both 
high and low redshifts that have metal absorption consistent with 
a cool, dusty environment. Sub-DLAs tend to have both higher 
metallicities and larger velocity widths than DLAs, and H2 is more 
likely to be found in DLAs with both these characteristics ( |Noter-| 
Idaeme et al.|2008l l. 

5.5 Evolution in 

We plot the /hj values as a function of cosmic time in Figure [To] 
There is no evidence for evolution in /hj , though more measure- 
was at a higher redshift, it may have been heavily blended with Lyo forest 
absorption. 
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Figure 10. The molecular mass fraction as a function of lookback time for 
sites where H2 absorption has been detected. The Local Group point is a 
median for all values from the plane of the Milky Way ( Savage et al. 1977|, 
the LMC and SMC l .Welty et al.|2012 l. The eiTor bars show the 10th and 
90th percentile level. There is no evidence for evolution in /h, , but the 
number of measurements at 2 > is small. 

ments are needed, particularly at intermediate redshifts, given the 
large scatter in seen along both local sightlines in the Milky 
Way halo and in higher z DLAs. 



6 SUMMARY 

We have analysed a sub-damped Lya system with A^hi^ 
■|^gi9.5±o.2 ^( 2 = 0.56 that shows associated molecular hy- 

drogen absorption in the Lyman and Werner bands. Using veloc- 
ity components determined from a high resolution spectrum cover- 
ing metal transitions falling in the optical, we fit a two-component 
model to the H2 absorption and find a lower limit to the molecu- 
lar fraction of logio /h2 ~ —1.93 ± 0.36, and a lower limit in- 
dependent of the assumed velocity structure of /h2 > 10"'^. This 
is higher than other sightlines with similar Nm where H2 has been 
measured in the Milky Way halo. We find a metallicity for the cloud 
logio Z = -0.86 ± 0.24, or 0.14l° J° solar. The dust-to-gas ra- 
tio relative to the solar neighbourhood is logio = —^-^-o't' or 
- = 0.04lHt 

We modeled the observed transitions using CLOUDY and were 
unable to find a single solution that can simultaneously reproduce 
all the observed transitions. However, a model for the absorber of a 
10* K cloud illuminated by a radiation field dominated by the UV 
background can broadly reproduce all the observed column densi- 
ties apart from those of H2, C I, and O VI. We conclude that there 
are three phases in the absorber; a T ~ 100 K phase where the C I 
and H2 arises, a T ~ 10* K phase where the low-ionisation metal 
absorption occurs, and a hotter, collisionally ionised phase associ- 
ated with O VI. 

Using simple models of H2 formation-dissociation equilib- 
rium, we calculate densities for the H2 absorbing region from 
~ 1 - 4cm"^ to ~ 70 - 480 cm , depending on the incident 
strength of the radiation field. The lower density range corresponds 
to cloud thicknesses of ~ 3 — 10 pc, the high density range to 
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~ 0.002 — 0.15 pc. Given the A^m, the presence of a three phase 
medium, the molecular fraction, metallicity and two galaxy candi- 
dates near the QSO sightline with impact parameters of ~ 10 kpc, 
we conclude this system may be a tidally stripped feature similar to 
the Magellanic Stream. 

Finally, we remark that of the seven sub-DLAs observed 
at jz < 0.7 for which there is the possibility to detect A^H2 It 
lO^^'^cm"^, two H2 detections were found. A survey for H2 in 
low-redshift sub-damped systems could be a fruitful way to mea- 
sure the physical conditions giving rise to these absorbers. 
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lar to this analysis were written using the NumPy and SciPy pack- 
ages (http : / /www. scipY . org), and plots were produced us- 
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Valid 


Apply to 






Segment 


A range (A) 


A setting 


a 


/3 


FUVA 


1600 - 1800 


1627 


1.233 X 10"^ 


-2.115 


FUVB 


1425-1600 


1589 


-0.952 X 10-3 


1.468 



Table Al. Wavelength shifts required to bring COS G160M exposures with 
central wavelength settings 1589 and 1627 A on to a common wavelength 
scale. The shift AA that must be added to a wavelength A is given by AA = 
aX + j3. The FUVA shift is applied to the 1627 central wavelength setting, 
and the FUVB shift is applied to the 1589 setting. 



A (A) 


AA (A) 


A 


AA 


1425 


-0.1009 


1625 


-0.0062 


1450 


-0.0848 


1650 


-0.0033 


1475 


-0.0679 


1675 


-0.0011 


1500 


-0.0514 


1700 


0.0009 


1525 


-0.0365 


1725 


0.0029 


1550 


-0.0245 


1750 


0.0052 


1575 


-0.0160 


1775 


0.0080 


1600 


-0.0101 


1800 


0.0116 



Table A2. Wavelength shifts inferred from the H2 absorption as function 
of wavelength. AA should be added at each wavelength A to correct the 
wavelength scale. 
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Zwaan M. A., Prochaska J. X., 2006, ApJ, 643, 675 



APPENDIX A: CORRECTING FOR WAVELENGTH 
SHIFTS IN THE COS EXPOSURES. 

When combining the COS exposures, we found there were shifts of ^ 
20 km s~^ in the wavelength solution between exposures taken using dif- 
ferent central wavelength settings. Table [aT] gives the shifts that must be 
applied to bring our G160M exposures with central w avele ngth settings 
1589 and 1627 A to a common wavelength scale. Table [a2| gives the fur- 
ther ~ 10 km s~^ shifts that were required to give an internally consistent 
wavelength solution based on the expected positions of H2 absorption fea- 
tures. 



